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ABSTRACT 
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The development of vehicles powered by direct-current 
sources together with the development of the silicon-control- 
led rectifier has led to the use of squirrel-cage induction 
motors operating on non-sinusoidal, variable-frequency volt- 
age supplies. A digital computer simulation of the tran- 


f. be a ati 
D haters Pass 
2 , i pathy 


Silent and steady-state pertormantéé of a three-phase motor 





is derived and its use in predicting motor operation is 


illustrated. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 


(All rotor parameters referred to the stator) 


V = Instantaneous value of voltage, volts 

di = Instantaneous value of current, amneres 

Lyy = Stator phase self-inductance, henry (assumed 
identical for all three phases) 

Ee, = Rotor phase self-inductance, henry (assumed 


identical for all three phases) 


Los = Stator-rotor phase mutual inductance maximum, henry 
C = Rotor position angle, radians 

1 = Rotor phase resistance ohms 

R = Stator phase resistance, ohms 

T = Electromechanical torque developed, newton-meters 
Tu = Mechanical torque, newton-meters 

po = Rotor angular velocity, radians/second 

J = Inertia of system referred to motor, kilogram meters 
Z = Friction and windage coefficient 

My = Magnetic flux linkages, weber 

p = d/dt 

A,B,C = Stator phases 

a,b,c = Rotor phases 

Ic = Time, seconds 

te = Time increment, seconds 


i iN TROBUEELON 

Engineers have long sought a method of combining the 
simplicity, compactness, low cost, ruggedness and low main- 
tenance characteristics of the squirrel-cage rotor type of 
induction motor with the speed control characteristics of 
direct-current motors, involving commutators and brushes, or 
of wound-rotor induction motors also using brushes. 

In recent years two developments have given importance 
and feasibility to this project. 

1) The present and projected develooment of space ve- 
hicles, deep-submergence-oceanographic vessels, and nuclear 
reactor coolant pumps has demonstrated and intensified the 
requirement for reliable, variable-speed drives which require 
no maintenance, not even brush replacement, during their ex- 
pected lifetime, and which are very reliable. Preferablv 
these drives should be able to operate from a direct-current 
source such as a fuel cell or battery pack. 

2) The development of the silicon-controlled rectifier 
and its use in bridge-inverters controlled by solid-state 
logic devices provides a reliable, practical method of ob- 
taining a variable-frequency three-phase power sunply from 
a direct-current source. These inverters commonly contain 
no moving parts, have very short switching times and anpear 
to be very reliable. 

Many descriptions of solid-state inverters have been 


[9] 


published in recent years and it will be assumed that 
voltage wave shapes of the following type are readily avail- 


able. 





Power Conditioner Voltage Output 


Figure 1-1 


70 ASSUMPTIONS 

1) The wave shapes of all the three phase voltages are 
cyclic and have the same shape. 

2) Iron losses are assumed to be small and are not 
accounted for. Saturation of the magnetic path is neglected. 

3) Deep-bar effects are not considered. 

4) Since the source voltage is applied to one phase in 
series with two phases in parallel (Y connection), and the 
windings are assumed to be symmetrical, the phase voltage 
division is one-third and two-thirds of the source voltage. 

5) The rotor is cylindrical as is the inside of the 
stator. The air gap is therefore uniform. 

6) The rotor is of the symmetrical, three-phase, coil- 
wound type. This is done to facilitate analysis although 
the basic phenomena occurring in a squirrel-cage rotor are 


the same. 
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ce PRINCIPLES OF OPERATION 





Single-Phase Motor Model 


Figure 3-1 


If at time zero a positive voltage v is applied to 


stator phase A as shown, current ly will flow as indicated 


creating flux which will link the turns of phase a. As ly 


increases an induced current will be generated in rotor 


phase a which will oppose the time rate of change of la: 


Fleming's right-hand rule may be used to establish the po- 
larity of this induced current as negative, i.e., in at the 
bottom and out at the top. By Lenz's law a force or torque 
is created which acts to turn the rotor in a counterclock- 


wise direction. Consider the sequence of events: Vv, Creates 


an exponentially rising current 1 1, produces flux linkage 


yey 


er the time rate of change of which generates voltage Va! 
and in turn v creates current i. 
It is obvious that if periodic wave forms are consid- 


ered 1. must lag i Also, since ve is created by the time 


ae 


rate of change of i,, when in begins to decrease in value 


A 


oy 


the polarity of induced voltage Ve must reverse. Therefore, 
at this point it may be assumed that the wave form of i, 
must lag the wave form of in: 
If now the rotor is considered free to rotate and it 

is assumed that counterclockwise torque is desired, the 
following relationships must be considered. Let 9, rotor 
angle, be zero when coil axes A and a are colinear and let 
counterclockwise be considered positive. Then the mutual 


inductance will equal Lx cos © Bbywdefinition. Torque, T; 


will then be equal to ig’ (-i,) (Se Lox cos 90] which equals 


*(-i_)-1b,, sin 6]. 
For positive torque with 0 greater than zero and less 
than II radians, the product in i, must be negative. For 0 


between II and 2II radians in i, must be positive. 

Another effect which must be considered is the rotor 
phase belt "carrying" its current around with it. In the 
numerical methods solution to follow it will be assumed that 
electrical transients are fast enough compared to mechanical 
transients so that 9 may be considered constant for the very 
short At periods involved. The foregoing does not mean that 
"speed voltages" or counter-electromotive force is neglected, 


but only that pO and mutual inductance are held constant for 


the At periods. 
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4, ANALYSIS OF THREE-PHASE MOTOR OPERATION 


Three-Phase Motor Model 


Figure 4-1 


If each coil is considered to have an axis along which 
a composite flux is created the total flux linkages along 
each axis may be computed as follows: 


Ye = ili + ipl. cos (2]1/3) + igh. cos (-2]1/3) + aay 


cos (0) + ipl cos (0+2II/3) + 1 dx cos (0-2]1/3) 


¥ 5 = inl, + om cos (2]1/3) + Tux cos (-2]1/3) + er 
cos(9-2]1/3) + ido. cos (@) + i Lay cos (9+2]1/3) 

Yo = ich, + ae cos (2]1/3) + Lal x cos (-2]1/3) + i vox 
cos (0421/3) + inlay cos(Q0-211/3) + i dy cos (0) 

= se + a eee cos (2711/3) + vax cos (-2]1/3 + 1 Loe 


cos (-0) + ce cos (2]1/3-0) + ighax cos (-2]1/3-0) 


b ipbys + See cos (-2][/3) + thx cos (2]1/3) + ae 


cos (-2]1/3-0) + ipbiy cos (-0) + ce eos (217 3-0) 
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ae = Le + tlle Gos (26). + x cos (-2]1/3) + Ls 


cos (2]1/3-0) + inl oy cos (-2]1/3-0) + LoL Gos (—G) 


The voltage equations for the physical machine are de- 
fined to be: 


Ve= IR + py 


Wee Le Rete od 


a A 11 Pta + Ls Gosent2il/ 3) pi, 


+ L COS (= 21773) pi, + L 


yx cos (0) pi, - i L Sg 


yx a yx 


(0) pO + tee cos (0+2I1/3) pi, - We Sin (0+2I/3) pQ 


7 Lo cos (0-2I1/3) pi, - ous Sin (90-211/3) pe 


Vp = ipR + b))pi, + ae cos (2I1/3) pig + ee cos (-2II/3) 


pi, + L cos (0-2II/3) Pileeac We Stn O=211/ 3) mpo 


Vx 


+ L COs eG) pi, - 1,4 


ve sin (0) pO + L cos 


yx yx 


(0+211/3) Pies Lhe Sin (0+2I1/3) pod 


Vo = igR + Ly, Pig + Lay cos (2II/3) pi, + Lao cos 


1 
ra 
ey 
CP) 
P- 
re | 


(-211/3) pint & cos (0+2I1/3) pi, 


yx 


(O0+2I/3) po + Ly cos (0-2]I/3) i Sale I sin 


COS ZIA 3 )aeeisy 4 Loose COs- (a) pi, - ly Hos sin (0) pod 


YS he + Loy pi. + ae Gos "(211 /3) Pip + te cos (-2II/3) 


SOs=( ©) pi, + are Sin (-0) pd 


+L cos (2II1/3-0) pi, + beth se Sin (2I1/3-0) pd 


+ L cos (-2]1/3-0) Pic + Toler Sin (-211/3-0) pd 


14 


ade 


ia > xx Pip * Lx 


Yb 
pl, + Lox 


ae Li 


yx os (-6) pi, + 


(211/3-0) pi, + Lobox 


er + Ln 
S mite Lyx Pr Lox 


Ply + box 


+ 1 COS 
yx 


(-O) pig + 


cos (-2II/3) 


sin 
ele see alilyaow) 
cos (2I1/3-0) pi, + tb 


(-21/3-0) pi, + igh 


pi. + Lx Cosco) 


cos (-211/3-0) pi, + Lalas Sin (-2[1/3-0) pd 


tts Sin (-0) pO + L cos 


Vx 


(211/3-0) pd 


pi, + Lx ees (-2]173) 


sin. (273-0) “pe 


Buy x sin (-211/3-90) p0 


Tone Sin (-0) p0 


If, for simplicity of notation, terms are redefined as 


follows; be= L cos 
yx 

E= L cos 
yx 

F = Lx cos 

Gs Lox cos 

He — ie cos 
yx 

ee =eb cos 
yx 

J = L cos 
yx 

re | cos 
yx 

a = & sin 
yx 

w= L sin 
yx 

Ee— | Sin 
yx 

QO = Lox sin 

oe 6 sin 
yx 

T= i sin 
yx 


ANS) 
(-211/3) 
(211/3-0) 
(211 /3+0) 
(-211/3=0) 
(-211/3+0) 
(0) 

(=@) 

(0) po 
(=O) Bo 
(211/3+0) pd 
(211/3-0) po 
(-211/3+0) pO 
(-21/3-0) po 


ike 


The equations in matrix form are: 





Li D E J G I ply 

E Loy D it J G Pl, 

D 1s Liy G 43 J Pla 

+ 

K F H L D E pl 
xX a 

H K F E ae D Ply 

| 

F H K D i bee Pl, 


It may be well to note at this point that no assump- 
tions or specifications have been made as to voltage wave- 
form or frequency. It has been assumed that the motor wind- 
ings are symmetrical, an assumption which simplified the 
notation and in addition allowed the use of the original 
voltage distribution assumption for a three-wire system. 

The modification required to consider a four-wire unbalanced 


system is obvious and will not be pursued. 
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The system of nonlinear differential equations obtained 
may be solved only if O and pO may be considered constant. 
Since the transient solution is desired this can only be jus- 
EEEtede ton very Shert time periods. 

The application of numerical methods and the use of a 
computer may be expected to produce an inexact but useful 
solution to the problem. 

Let the previously derived system of matrix equations 
be expressed by: 

(v] = [H] (iJ + ([K] [pi] 
then: 
[pil = (K]-"[v) - (k]7> tH) (4) 

Torgue generated in the constant air-gap machine is 
equal to the summation of the product of stator phase cur- 
rents, rotor phase currents, and the space rate of change in 


their respective mutual inductances. 


Therefore: 
Ta ~i,i bay sin © - iniploy Sin (0+2I[/3) - A ely. 
sin (0-211/3) - she ee sin 0 - toby Sin (0+2II/3) 
- a au Sin (0-2I[/3) - Loi chy sin 0 - Lotalys 


sin (0+211/3) - doit yy Sin (0-21/3) 


The mechanical equation for torque: 


T = T., = J(p(p0)) + Zpod 


may be solved for p(p0): 


p(po) = (T-T_-Zp0) /J 


Py 


If load torque is assumed to be known, predictor- 


corrector methods may be used as follows: 


PREDICTOR 
(pi), = £(0,p9,i),, 
p(po), = £(i,9),, 
14] i (At) (pi), 
CORRECTOR 
(pi), = £(0,p0,1) 44 
Ioel = I, + dt (pirPI) 72 
Gheek 
et ae 
EEEOr = = 
—n+1l 


If the absolute value of error is found to be unaccept- 


able, set I, = I and recompute i. If the error val- 


+1 =—nt+1 =pll 
ues are acceptable, torques, rotor angles and angular veloc- 
ity may be determined preparatory to the next step. Adjust- 
ing At after each cycle of computations will permit saving 
of computer time while maintaining the desired degree of 
accuracy. 

Note that all motor parameters used in determining ex- 
pected transients or steady-state response may be determined 
by relatively simple laboratory tests of a given machine. !7! 
Stator and rotor copper (te) loss calculations may be in- 
Corpenbaued Ineo the program. 

TE us theoretically osemple to account for hysteresis 
and eddy-current losses by detailed analysis of flux changes 


and the magnetic path. However, the parameters required for 
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these calculations are not normally available to the engineer 
for am "off the shelf" motor. Since the induction motor 
circuits themselves are heavily inductive it may optimis- 
tically be assumed that the circuits act as a low-pass filter 
and typical loss figures for the same class of motor overa- 
ting on a three-phase balanced sinusoidal voltage supnnvly may 
be used. 
Su APPLICATION 

For purposes of illustration a 20-horsepower, 440-volt, 
three-phase, 60-Hz., 3470-rpm, two-pole, squirrel-cage in- 
duction motor is selected. The stator windings are WYE con- 
nected and the phase model parameters, as determined by 


[5] 


blocked-rotor and no-load tests are as follows: 


R = 0.30 ohms 


r = 0.36 ohms 
L, = L, = 0.00167 henrys (leakage inductances) 
Le model = 0.0723 henrys (magnetizing inductances) 


The model magnetizing inductance being three-halves the 


actual value, the inductance values used are: 


Lii= Ly + Lo (273). = 0a? 

Liter Loi = 0.04987 

L = L model (2/3) = 0.0482 
yx yx 


The mechanical load is arbitrarily assumed to be a pro- 
peller with load torque equal to 0.000315 times the square 
of the angular velocity in radians per second. This coeffi- 


cient is selected to require 19.5 horsepower delivered to 


ro 


the load at rated speed and the friction coefficient (0 .@iie2) 
accounts for the remaining one-half horsepower at rated con- 
aiewonsi. 

The program following, Appendix I, was successfully run 
on an IBM 360 computer. Limiting the storage requirements 
to that required for 200 steps at a time allowed double- 
precision computation using less than 100,000 bytes of core. 

IBM supplied the subroutine for matrix inversion (GAUSS3) . 
Provision has been made to avoid multiplication or division 
by near-zero or zero values. 

SINUSOIDAL VOLTAGE SUPPLY 

In order to obtain a standard of comparison, a three- 
phase, 60-Hz, sinusiodal voltage supply was assumed for simu- 
lation purposes. Figures 6-1 through 6-6 display the com- 
puted current wave-forms and generated torque. The areas of 
special interest are considered to be motor starting tran- 
Sients and steady-state performance. 

NON-SINUSOIDAL VOLTAGE SUPPLY 

A 60-Hz, pulsed voltage supply is assumed as an output 
from the power conditioner which, in accordance with assumpn- 
tion number four, is resolved in the ratio one-third and 
two-thirds across the stator phases. Comvouted phase voltage 
wave-forms, current wave-forms and generated torque values 
are plotted against time for starting and steady-state con- 
ditions in Figures 6-7 through 6-10. Voltage magnitudes are 
computed to have the same rms value as the rated sinusoidal 


voltage. 
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VARIABLE-FREQUENCY VOLTAGE SUPPLY 

Assumptions are the same as above except that tachom- 
eter feedback is assumed to control the frequency and thus 
the pulse widths of the stator voltage supply. If the cir- 
cuit model Thévenin equivalent is solved as a function of 
phase-voltage frequency and rotor angular velocity for maxi- 
mum power transfer to the load, and thus maximum generated 
torque, and if stator-phase resistance is neglected, the op- 
timum voltage frequency is found to be equal to rotor angu- 
lar velocity plus a constant. The results of this simulation 
are recorded in Figures 6-11 through 6-13. 

While very large starting torques are obtained, the 
currents generated are excessive and any motor not especially 
designed would probably be destroyed by a short neriod of 
Operation involving a series of mechanical transients. The 
current requirements and heat dissipation requirements for 
the power conditioner would also far exceed those required 
for steady-state operation. 

Consider the steady-state relations that magnetic flux 
linkage per pole is proportional to the ratio of rms voltage 
to frequency and that torque generated is proportional to 


[5] 


the space rate of change of magnetic flux linkages. Con- 
sider also that the variable-frequency pulse width described 
above is one-sixth of the stator-voltage cycle time. The 

rms voltage for a given frequency may be controlled by varv- 


ing the pulse width up to a maximum of one-sixth of the cycle 


time. If the pulse width is made proportional to frequency, 


Zn 


the conditions quoted above are satisfied; therefore, maximum 
rated torque will be available during mechanical transients, 
and rms currents will be reduced. The results of this simula- 
tion, including the described voltage wave-form, are recorded 
as Figures 6-14 through 6-17. 
Or CONCLUSIONS 

Results displayed in Figures 6-1 through 6-13 compare 


[8] 


favorably with those obtained by other methods. Figure 
6-18 is a comparison of pO vs. time for a given motor opera- 
ting on different voltage wave-forms. 

From the data obtained the following conclusions may be 
drawn: 

The pulsed wave-form output obtainable from a power 

conditioner is useable for squirrel-cage motor 

operation. 

Variable speed control may be obtained by variation 

of power conditioner frequency. 

If the load requirements are such that few large 

speed changes are expected, such as normal opera- 

tion of a nuclear reactor coolant pump, the added 

complexity of the power conditioner required to con- 

trol pulse-width will probably not be justified. 

If large and frequent speed changes are expected, 

as would probably be the case in provelling an ocean- 

ographic vessel, the added complexity involved in 


pulse-width control is justified. 
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The computer program used to simulate the last wave-form 
described is presented as Appendix I. Computations using 60 
time increments per stator voltage cycle required anprox- 
imately 11 minutes of IBM 360 Computer time to simulate two 
and one-half seconds of motor transient time. 

Areas of further investigation should include considera- 
tion of the effects of magnetic saturation and its effect on 
current and generated torgue. Inclusion of power conditioner 
operation, copper losses, and vibration calculations in the 
simulation would improve the program's validity and thus its 


usefulness. 
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Stator Phase Voltages vs. Time. Constent- 
Frequency, Non-Sinusoidal Voltage Supply. 
Motor Undergoing Starting Transients. 
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APPENDIX I 
FORTRAN IV COMPUTER PROGRAM FOR SIMULATION OF THREE PHASE 
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